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ABSTRACT 

HH 212 is a nearby (460 pc) protostellar jet discovered in H2 powered by a Class source, IRAS 
05413-0104, in the L1630 cloud of Orion. It is highly collimated and symmetric with matched pairs of 
bow shocks on either side of the source. We have mapped it in 850 /im continuum, SiO ( J = 8 — 7), CO 
(J = 3-2), SO (Nj = 8 9 -7g), HCO+ (J = 4-3), andH 13 CO+ (J = 4-3) emission simultaneously at ~ 
1" resolution with the Submillimeter Array (SMA) . Thermal dust emission is seen in continuum around 
the source, mainly arising from an inner envelope (i.e., the inner part of a previously seen flattened 
envelope) and a possible disk. The inner envelope is also seen with rotation in CO, HCO + , and probably 
SO. Like H2 emission, CO and SiO emission are seen along the jet axis but extending closer to the source, 
tracing the bow shocks with a broad range of velocities and the continuous structures in between. SO 
emission is seen only around the source, forming a jetlike structure extending along the jet axis from the 
source, likely tracing the jet near the launching region. The jet is episodic and bending. It may also be 
slightly precessing as the jetlike SO structure shows a slight S-shaped symmetry about the source. A 
hint of jet rotation is also seen across the jet axis. Internal outflow shells are seen in CO and HCO + , 
associated with the bow shocks in the inner part of the jet. The bases of the HCO + shells are seen with 
a hint of rotation similar to that seen in the inner envelope, probably consisted mainly of the material 
■^j- ' extended from the inner envelope and even the possible disk. The bases of the outflow shells are also 

seen in H 13 CO + and even the continuum, probably tracing the dense material extended from around the 
same regions. Outflow shells are also seen in CO surrounding the jet extending out from the H2 nebulae 
seen around the source. 

Subject headings: stars: formation — ISM: individual: HH 212 — ISM: jets and outflows. 
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1. introduction jet w ith proper motions along the jet axis ([Claussen et alj 

On 1 Protostellar jets are seen associated with low-mass pro- ' ™ e & interacts with the ambient material, driv- 

1 ■ , , . , , c , r . ■ t -4- c m S a collimated CO outflow around it (Lee et al. 2000: 

^ tostars m the early stages of star formation. In spite of T , . „„„„ , „, a N „ v , , . . 

\J , . i . , ,. / , Lee et al. 2006, hereafter Paper I). Deep observations 111 

s_i ■ numerous studies, their physical properties (e.g., speed, TT . , . ' . „ — , — ; . K . „ 

, . j. „. ,. j 1 , \ H9 using the Infrared Spectrometer And Array Camera 

episodic nature, colhmation, and angular momentum) are , T „ . . J? „ , . ™^ ,. r . . ,,. T i, . 

still not woll understood (Hartiean et al 20001 Thcv are (ISAAC) of the ESO Very Large Telescope VLT show a 
Co SL111 not wen uiiticiSLOOQ inaingan ll ai. zuuui. _Liicy aie . r ltrr , , , r . , , , , 

, ,. j j. v 1 v j r 7; T. 1 j pair of diffuse nebulae near the bases of the iet, probably 

^ believed to be launched from accretion disks around the f ,, , „ ., „ ... 



protostars (|Shu et all F2OOOL iKonid & Pu dritz 2000), al- * racin S the outflo ^ ^ vit ^ walls illuminated by the bright 

; , ,/ , . , . bow shocks around the source McCaughrcan et al. 2002, 

lowing us to probe the accretion process, which remains . r , ,,„.,„.„- _,. . ' * — — ; r~~ . — , 

, , c i j • . , , ,. hereafter MZAML02 . I he inner part of the iet is also 

heretofore unresolved, as it requires us to observe directly . , , ; „.„ . . \tt, — -J — n Un»^i 

rrt ■ • 4.1 • 4. 4.4.1 att 1 mv. o u -it 4. a seen m shock excited SiO emission Chapman et al. 2002: 

L « , m the inner parts at the AU scale. Ihc Submillimeter Ar- . — ■ — rn u^^ a ittt - ; : — ; — n u„AA — , , . , , — r~ 

ray (SMAfl dHo et al.ll2004 . with the capability to probe ^ et . aL | l 20 ,^ \ Ta f™ et al ™ f . and ^J^J 

, j r - ^ , , , . 1 114.- CO emission (Paper 1). Since the iet is compact, dense, 

warm and dense molecular gas at high angular resolution. . , .\ — B — fr. . f. . ... ' .' 

1 j4- 4.j4.iu -i 4.- and warm, high-resolution observations m high-J transi- 

can be and has been used to study the physical properties , . .. Vn-^ 1 ,1 11, , 1 

c 4.1 . 4. . 4. j 4. -l ITT- — iilonnd RTi — — Tl tion lines 01 SiO, CO. and other shock tracers are presented 

ol the icts m great detail Hirano et al. 2006; Palau et al. . , ,. 

2006]) re s ^ u( ly its physical properties. 

The HH 212 jet, which is nearby (460 pc) in the L1630 
cloud of Orion and remarkably well-defined, is one of the 2. observations 

best candidates to investigate the physical properties of Observations toward the HH 212 jet were carried out 

the jets. It was discovered m shock excited H 2 emission with the gMA Qn DeC ember 2, 2005 in the compact config- 

dZmnecker, McCaughrean fc Rayner i| 1998| ), powered by a uration and on January 14; 2 006 in the extended configu- 

low-lummosity (~ 14 L ) Class source (or protostar), ration SiQ (J = 8 _ ?); co (J = 3 _ 2)j SQ ( Nj = gg _ 7g); 

IRAS 05413-0104. It is highly collimated and highly sym- HCO+ ( J = 4 — 3) and H 13 CO+ ( J = 4 — 3) lines were 

metric with matched pairs of bow shocks on either side observed simultaneously with 850 /an continuum using the 

of the source. Water masers are seen at the bases of the 345 GHz band receivers The rest frequency, upper energy 
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level, and critical density of these lines are given in Ta- 
ble [T] The receivers have two sidebands, lower and upper, 
covering the frequency range from 345.5 to 347.5 and from 
355.5 to 357.5 GHz, respectively. Combining the two side- 
bands results in a total bandwidth of 4 GHz centered at 
351.5 GHz (or A ~ 850 /im). Seven antennas were used in 
the array, giving baselines with projected lengths ranging 
from 8.5 to 180 m. Since the primary beam has a size of ~ 
35", two pointings with a separation of 18" were used to 
map the inner part of the jet within 25" from the source. 
For the correlator, 256 spectral channels were used for each 
104 MHz chunk, resulting in a velocity resolution of ~ 0.35 
km s _1 per channel. 

The visibility data were calibrated with the MIR pack- 
age, with Saturn, Uranus, and Quasar 3C454.3 as pass- 
band calibrators, Quasar J0530+135 as a gain calibrator, 
and Uranus as a flux calibrator. The calibrated visibil- 
ity data were imaged with the MIRIAD package. The 
dirty maps that were produced from the calibrated visi- 
bility data were CLEANed using the Steer clean method, 
producing the CLEAN component maps. The final maps 
were obtained by restoring the CLEAN component maps 
with a synthesized (Gaussian) beam fitted to the main 
lobe of the dirty beam. With natural weighting, the syn- 
thesized beam has a size of l"16x0"84 at a position angle 
(P. A.) of 38°. In some cases, however, in order to study 
the compact structures (i.e., the jet and the possible disk) 
in the system, the visibilities with uv radius shorter than 
40 m (or > 2"5) are excluded, resulting in a smaller syn- 
thesized beam with a size of 0"96x0"69 at P.A.= 38°. The 
rms noise level is ~ 0.16 Jy beam -1 in the channel maps 
and 5 mJy beam -1 in the continuum map. The veloci- 
ties of the channel maps are LSR. The absolute positional 
accuracy in our maps is expected to be one tenth of the 
synthesized beam, or ~ 0"1. 

3. RESULTS 

In the following, our results are presented with the H2 
image made with the VLT (MZAML02) for comparison. 
The H2 image, which shows a clear jet consisting of a 
chain of bow shocks and sinuous continuous structures 
in between in great detail, is used as a reference for the 
jet. Since our observations were carried out about 4 years 
later than the H2 image, the bow shocks, which were 
found to be moving at 100-200 km s" 1 or 0"045-0"09 
per year (MZAML02), might have moved down along the 
jet axis by 0"18— 0"36 from their positions in that im- 
age. This movement of the bow shocks, however, does 
not affect significantly our comparison and conclusions, 
considering that the angular resolution of our observa- 
tions is ~ 1" along the jet axis. The H2 jet is almost 
in the plane of the sky, with the blueshifted side to the 
north and the redshifted side to the sout h of the source 
(Zin necker. McCaughrean. k, Ravneril 19981 ). The systemic 
ve locity in this region is assumed to be 1.7 ±0.1 km s _1 , as 
in iPaper J . Throughout this paper, the velocity is relative 
to this systemic value. 

3.1. 850 fim Continuum Emission 

A continuum source is detected at 850 /xm in between 
the diffuse H2 nebulae with an emission peak at a(2ooo) = 
05 h 43 m 51?404 ± 0"1, 5 (20 oo) = -01°02'53;'10 ± 0"1, and 



a total integrated flux of 0.38±0.10 Jy (Fig. QJi). This 
peak position is within 0"1 from that found at A = 
3.5 cm with the VLA at an angular resolution of 0"3 
(|Galvan-Madrid e t al. 2004]), and is thus considered as the 
source position throughout this paper. The structure of 
the continuum source is not well resolved in the restored 
map. It is better seen in the CLEAN component map, 
which shows a faint flattened structure with a radius of ~ 
1" (460 AU) perpendicular to the jet axis and a bright un- 
resolved compact structure with a radius < 0"3 (138 AU) 
at the center. Faint emission is also seen extending to the 
northeast and southwest. The faint flattened structure, 
which is not resolved in the minor axis, may trace the in- 
ner envelope, i. e., the in ner part of the edge-on flattened 
envelope seen in lPaperll , while the compact structure may 
trace the envelope further in and a possible disk. The faint 
emission extending to the northeast and southwest, on the 
other hand, may trace the material extended from the in- 
ner envelope, as discussed later in £14.21 A compact source 
with an envelope is also seen in the amplitude versus uv 
distance plot. However, observations at higher angular 
resolution are really needed to confirm the structure seen 
in the CLEAN component map. The spectral energy dis- 
trib ution (SE P) of the continuum source (see Fig. [Tb and 
also lPaper J ) indicates that the continuum emission at 850 
/jm is mainly thermal dust emission. Assuming a constant 
temperature, a mass o pacity n v = 0.1(^/10 12 Hz)' 3 cm 2 g" 1 
(|Beckwith et al.lll990h , and a source size of 1 arcsec for 
the dust, the SED can be fitted with = 1 and a temper- 
ature of ~ 48 K. Thus, t he (gas + dust) mass is estimated 
to be ~ 0.08 M e , as in IPaper 1 With 0=1, the bright 
compact structure at the center may inde ed harbor a disk 
around the source (|j0rgensen et al.ll2006f ). 

3.2. Jet Axis and Bending 

The northern component and southern component of 
the H2 jet are not exactly antiparallel (Figs. and [2b). 
Their inner parts are misaligned by ~ 2°, with their axes 
found to have a P.A. of 21.5°±0.5° and 203.5°±0.5°, re- 
spectively, by connecting the source to the inner H2 bow 
shocks out to bow shocks NK7 and SK7. Since their orig- 
inal paths of motion are likely to be antiparallel, this mis- 
alignment suggests a presence of a jet bending (see H4.1.6I 
for possible mechanisms). Assuming both bent by the 
same degree, their inner parts are both bent by ~ 1° to 
the west with the original axis having a P.A. of 22.5°±1°. 
Note that the bending actually continues further out and 
beyond bow shocks NB1/2 and SB1/2. 

The outer parts of the jet may have different original 
axis and bending because the tips of bow shocks NB3 and 
SB4 are on the opposite sides of the axes of the inner parts. 
The original axis, which can be estimated by connecting 
the tips of bow shocks NB3 and SB4, is found to have 
a P.A. greater than that of the inner parts, suggesting 
that the jet axis might have rotated clockwise more than 
a few thousands years ago. However, without knowing the 
source position at the time when bow shocks NB3 and SB4 
were formed, we are not able to determine the bending. 

3.3. SiO, CO, and SO emission 
3.3.1. Morphologies 
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SiO emission is detected mainly along the jet axis 
(Fig. [2]:). About 70% of the single-dish flux is recov- 
ered from our observations, comparing to the Atacama 
Submillimeter Teles cope Experiment (ASTE) observations 
(jTakami et al.f 2006). ft is detected toward the bow shocks 
and the continuous structures seen in H2. Bow shocks 
NK2, NK4, SK2, and SK4, which are bright in H 2 , are 
also bright in SiO. However, bow shocks NK1 and SKI, 
which are bright and prominent in H2 , are faint in SiO and 
only detected on their eastern sides. SiO emission is also 
detected at ~ 1" north and south of the source, forming 
knotty structures labeled knots SN and SS, respectively 
(see Fig. for the closeup). The knots are spatially un- 
resolved with a diameter (transverse size) of less than 1" 
(460 AU) . A faint elongated structure is seen extending to 
the southwest from knot SS, fitting in between the wings of 
bow shock SKI, possibly tracing the "intrinsic" jet itself, 
but not the shock. 

CO emission is detected not only along the jet axis but 
also toward outflow shells (Fig. [4ji). In the following, two 
velocity ranges, high and low, are selected to show these 
two components separately. At high velocity (from —18.4 
to —7.1 km s" 1 and from 3.4 to 13.3 km s _1 ), CO emission 
is seen mainly along the jet axis toward the bow shocks 
and the continuous structures seen in H2 , with some coin- 
cident with the SiO emission (compare Figs. |2j: and[2jl). 
CO emission is also detected at ~ 1" away from the source 
(see Fig. [3b for the closeup), associated with knots SN and 
SS seen in SiO. At low velocity (from —4.0 to 3.1 km s _1 ), 
limb-brightened outflow shells (labeled A and B) are seen 
extending to the north and south, respectively, from the 
H2 nebulae (Fig. |4Jd). They seem to connect to the H2 bow 
shocks SB1/2, NK7, and NB1/2 and some further down 
the jet axis, and are thus probably internal outflow shells 
driven by them. Throughout this paper, outflow shells are 
called internal if not driven by the leading bow shocks at 
the heads of the jet. Near the source, a limb-brightened 
internal outflow shell is seen in the south connecting to 
bow shock SKI (Fig. 0J:). Internal outflow shells are also 
seen in the north but associated with more than one bow 
shock. The emission extending along the jet axis to bow 
shock NK1 is unresolved but may trace an internal outflow 
shell as a counterpart of that extending to bow shock SKI. 

SO emission, on the other hand, is detected mainly 
around the source (Fig. [3};). Jetlike emission is seen along 
the jet axis extending from the source to the north and 
south to knots SN and SS, respectively, likely tracing the 
jet near the launching region. The northern component is 
shifted slightly to the east and the southern component is 
shifted slightly to the west of the jet axis, showing a slight 
S-shaped symmetry about the source. 

3.3.2. Kinematics 

The kinematics of the bow shocks and the continuous 
structures can be studied with the position-velocity (PV) 
diagrams of the SiO and CO emission cut along the jet axis 
(Fig. [5]). In the PV diagrams, the CO emission around 
the systemic velocity, which merges with that of the ambi- 
ent cloud, and the CO emission from 7 to 9 km s _1 , which 
merges with that of the foreground ambient cloud with the 
same vel ocity ran ge, are resolved out from our observations 
(see also iPaper t ) . The internal outflow shells (labeled S) 



near the source are seen in CO with the velocity magni- 
tude increasing with the distance from the source toward 
bow shocks NK1, NK2, and SKI on both the redshifted 
and blueshifted sides, indicating that the shells around the 
bow shocks are expanding mainly transversely perpendic- 
ular to the jet axis. This is expected f or outflow shell s 
driven by jet-driven bow shocks (see, e.g.. lLee et al.ll200~lh . 
The bases of the shells, on the other hand, may also have 
velocity component parallel to the jet axis. 

Like the H2 emission, the SiO and CO emission toward 
the bow shocks and the continuous structures are mainly 
blueshifted in the north and redshifted in the south. They 
are, however, more blueshifted in the north than redshifted 
in the south, with a mean (representative) velocity of ~ 
— 11 km s _1 in the north while 7 km s -1 in the south, sim - 
ilar to that seen in the H2 emission (jTakami et all 12006(1 . 
The reason for this asymmetry in the velocity is unclear. 
The SiO emission is seen with a range of velocities toward 
the bow shocks: bow shocks NK2, NK4, SK2, and SK4 are 
bright and seen with a velocity range of ~ 15 km s _1 , while 
other bow shocks are faint and seen with a narrower veloc- 
ity range. The SiO emission is not well resolved, probably 
some from the bow tips with a broad range of velocities 
and some from the bow wings with the velocity magni- 
tude decreasing rapidly away from the bow tips. The CO 
emission is also seen with a range of velocities toward the 
bow shocks, some with similar and some with narrower 
velocity range than that of the SiO emission. Some CO 
emission is probably from the bow tips and some from the 
bow wings, with some (e.g., bow shocks NK2 and SK2) co- 
incident with the SiO emission. Note that a broad range 
of velocities is also seen ahead of bow shock NK4 in CO, 
suggesting that a shock is formed there. 

The continuous structures between the bow shocks may 
trace the "intrinsic" jet itself and thus can be used to study 
the jet kinematics. In the south, the continuous structure 
is seen between bow shocks SKI and SK2 in CO and SiO. 
The SiO and CO emission there together may show that 
its velocity increases with the distance from bow shocks 
SKI to SK2. In the north, the continuous structures are 
seen between bow shocks NK1 and NK2 and between bow 
shocks NK2 and NK4 in CO. Their velocity also seems to 
increase with the distance from one bow shock to the next. 

The kinematics of knots SS and SN around the source 
can be studied with the PV diagrams of the SiO, CO, 
and SO emission (Fig. In SiO, knots SN and SS 

are also seen with a broad range of velocities, but with a 
lower mean velocity than that seen toward the bow shocks 
further away. The kinematics of the CO emission is un- 
clear, probably with a velocity magnitude increasing with 
the distance from the source toward the knots on both 
blueshifted and redshifted sides. The kinematics of the 
SO emission is also unclear, probably partly similar to 
that of the SiO emission and partly similar to that of the 
CO emission. Observations at higher resolution are really 
needed to resolve their structures and study the detailed 
relationship among them. 

3.4. HCO+ and H 13 CO+ emission 

3.4.1. Morphologies 

HCO + emission is seen extending to the north from the 
source surrounding the CO emission that extends to bow 
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shock NK1 (Fig. [7^,). The emission is likely from a limb- 
brightened shell as seen in the blueshifted emission, but 
with the east side brighter than the west side (Fig. [7b>). It 
extends to the faint bow shock NF and then to the wings 
of bow shock NK1. Its base is also seen in continuum co- 
incident with its northeast extension (Fig. [7b). In the 
south, a similar shell is expected to be seen in HCO + , but 
it is faint and only seen at the base surrounding the CO 
shell that extends to bow shock SKI. As in CO, the shell 
is brighter on the west side. Note that the emission (la- 
beled NW in Fig. [7J:) extending to the northwest from the 
source is almost perpendicular to the jet axis coincident 
with the continuum emission and thus may arise from the 
inner envelope. 

A shell is also seen in H 13 CO + extending to the north 
from the source, surrounding the east side of the north 
base of the HCO + shell (Fig. [7Ji) . The emission seen 
across the source may arise from both the inner envelope 
and the shell, with the emission peak northwest of the 
source. 

3.4.2. Kinematics 

As mentioned, the HCO + emission around the source 
perpendicular to the jet axis may arise from the inner en- 
velope. A cut across the source perpendicular to the jet 
axis shows that the emission is seen across the source with 
the redshifted and blueshifted peaks on either side of the 
source (Fig[5^), similar to that seen in the C 18 and 13 CO 
emission (|Paper J ) but closer to the source tracing the in- 
ner envelope. The blueshifted emission has a peak at ~ 
(— 1 km s , — 0"05) in the west, while the redshifted emis- 
sion has a peak at ~ (1 km s~ 4 , 0"05) in the east. This 
PV structure suggests that the inner envelope is not only 
rotating around the source but also infalling toward the 
source with the redshifted emission from t he nearsi de and 
the blueshifted emission from the farside (|Paper J ). The 
emission is optically thick, with the absence of the emission 
around the systemic velocity. The blueshifted emission 
is much brighter than the redshifted emission, co nsistent 
with a presence of an infall motion (|Evansl[l999l ). Note 
that the redshifted emission in the west at ~ (0.7 km s _1 , 
—2") may arise from the shell. On the other hand, the 
origin of the H 13 CO + emission is unclear and may arise 
from both the inner envelope and the shell. 

The kinematics of the shells can be studied with the cuts 
across the shells centered at the jet axis. A cut across the 
shell at 2"5 north of the source shows two opposite HCO + 
crescents at ~ ± 1 km s" 1 and one H 13 CO + structure 
around the systemic velocity (Fig. [5b) . The two crescents 
may form the two ends of an elliptical structure connected 
by the H 13 CO + structure, suggesting that the shell there 
is expanding at ~ 1 km s~ 4 . The HCO + shell is probably 
optically thick, with the absence of the emission around 
the systemic velocity. The redshifted emission is much 
brighter than the blueshifted emission, consistent with a 
presence of an outflow motion. Cuts across the shell at 
l'/25 north (Fig. [8b) and 1"25 south (Fig. [SU) of the 
source show a hint of rotation with the redshifted emis- 
sion to the east and the blueshifted emission to the west, 
similar to that seen in the inner envelope, suggesting that 
the shell is probably the material extended from the inner 
envelope. Note that, however, the center of symmetry is 



shifted by ~ 0''5 to the east in the north, while 0"5 to the 
west in the south from the jet axis. 

3.4.3. Temperature, Column density, and Density 

In the following, we derive the excitation temperature, 
column density, and density, assuming optically thin emis- 
sion and local thermal equilibrium. Note that due to the 
absorption of the ambient cloud, that part of the emission 
is optically thick, and that part of the emission is resolved 
out by the interferometer, the values of the column den- 
sity and density presented here are lower limits of the true 
values. 

The excitation temperature of the CO emission can be 
derived from the line ratio of CO J = 3 — 2/ J = 2 — 1 along 
the jet ax is (Fig. H?h.), using the CO J = 2 — 1 observa- 
tions from lPaper 1 Note that, however, since the emission 
around 7 km s _1 is resolved out more in CO J = 2 — 1, 
the line ratio there is not properly calculated and appears 
to be larger than the maximum value of 2.25. In addition, 
since some of the emission, mostly at high velocity, is too 
weak to be detected in CO J = 2 — 1, the line ratio there 
also appears to be larger than the maximum value of 2.25. 
At a given projected distance from the source, the line ra- 
tio increases from low velocity to high velocity, suggesting 
that CO J — 3 — 2 traces high velocity (and thus jet) better 
than CO J = 2 — 1. The temperature of the shells is found 
to be ~ 20 K (Fig. [5b>). The temperature of the high- 
velocity emission in the jet is higher and assumed to be 50 
K, with a peak brightness temperature found to be ~ 35 K 
in knot SS. The H2 column density can be derived assum - 
ing a CO abundance of 8.5xl0~ 5 (jFrerking et al.lll982D . 
In the jet, it is found to be ~ 0.4xl0 21 cm~ 2 in the con- 
tinuous structures and ~ 1.2xl0 21 cm~ 2 in knots SS and 
SN. Therefore, with a diameter of ~ 1", the continuous 
structures has a density of 5.8xl0 4 cm -3 , and knots SS 
and SN have a density of 1.7 xl0 5 cm" 3 , similar to tha t 
found in the CO jet of HH 211 (|Gueth fc Guillotead[l999h . 
On the other hand, the shells are found to have a column 
density of (2 — 7) x 10 20 cm~ 2 and thus a density of ~ 
(0.3 - 1.0) x 10 5 cur 3 , with a shell thickness of ~ 1". 

The excitation temper ature of the SiO e mission has been 
found to be 50-150 K (jGibb et al.ll2004h and is thus as- 
sumed to be 100 K here. This temperature, derived at 
much lower angular resolution, can be considered as the 
lower limit of the true value. With this temperature, the 
SiO column density is found to be 1.5 — 2.5 x 10 14 cm~ 2 
toward the bow shocks and k nots. The SiO ab undance 
has been found to be 5 x 10~ 8 (|Gibb et al.l l2004). With a 
size < 1", the density is > (4 - 7) x 10 5 cm -3 . Note that, 
since the temperature of the SiO e mission could be hi gher, 
i.e., 300-500 K as seen in HH 211 (jHirano et al.ll2006f) . the 
density could be a factor of a few higher. 

The excitation temperature of the SO emission can be 
derived from the line ratio of SO Nj — 89 — 1%/Nj = 
pa — 4,5 , using the SO Nj — 5q — 4 5 observations from 
iPaper J . With a line ratio of ~ 1.5, the excitation temper- 
ature is found to be ~ 100 K, similar to that found in the 
CepA-East outflows (jCodella et alJ feOOS). The SO column 
density is found to be ~ 8xl0 14 cm~ 2 . In shock model, 
the SO ab undance is predicted to be 3 x 10~ 8 — 6 x 10 -7 
at 100 K (|Viti. Nataraian. fc Willia ms 2002). Thus, with 
a size of ~ 1", the density is 2 x 10 5 - 4 x 10 6 cm -3 . 
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The excitation temperature of the HCO + emission is 
close to the dust temperature, with a peak brightness tem- 
perature found to be ~ 41 K toward the source. Assuming 
a temperature of 50 K, the HCO + column density is found 
to be ~ 3.2 xlO 13 cm -2 toward the source and 1.8xl0 13 
cm -2 toward the shells. The abundance of H CO + is uncer- 
tain. It is ~ 2 x 10 ~ 9 in molecular cloud cores (jGirart et alj 
2000j) but greatly enhanced in outflow interacting regions 
by a factor of 20 with a value of 4x 10~ 8 (Hogerheiide et al.l 



Il998t IViti. Nataraian. fc Williamsll2002f) . Therefore, with 
a thickness of ~ 1", the density is 1.2 x 10 5 - 2.4 x 10 6 
cm -3 toward the source and 6.5 x 10 4 — 1.3 x 10 6 cm" 3 
toward the shells. 

The excitation temperature of the H 13 CO + emission 
is unknown and assumed to be the same as that of the 
HCO + emission. The H 13 CO + column density is found to 
be ~ 6 x 10 12 cm -2 . The H 13 CO + abundance has been 
found to be ~ IP -10 . m molecular cloud cores (see, e.g., 
iTakakuwa et al. 2003), but is unknown in outflow inter- 
acting regions and could be enhanced by a factor of 20, as 
the HCO+ abundance. Thus, the density is 4x 10 5 -8x 10 6 
cm -3 , with a thickness of ~ 1". 

4. DISCUSSION 

4.1. Molecular Jet 

4.1.1. Shocks: SiO, CO, and H 2 

SiO is seen trac ing the shocks in the jet, as in HH 211 
(Pal au et al.ll2006f h The emission is seen around the bow 
shocks with a range of velocities, as predicted in the jet- 
driven bow shock models (see, e.g., iLee et al.ll200l . Knots 
SN and SS are likely unresolved bow shocks, associated 
with a similar range of velocities. H2 emission is also seen 
associated with them in the Spitzer IRAC 4 (8 /im) and 
IRAC 3 (5.8 /im) images. No H2 emission is seen asso- 
ciated with them in the VLT image at 2.12 /im is likely 
because of the severe dust extinction toward them at that 
wavelength. The spectra toward bow shocks NK4, SK2 
and SK4 and knot SS are asymmetric: steep toward the 
highest velocities with a wing pointing toward the lowest 
velocities (Fig. [TO]) , also as predicted in the shock model 
in which the S iO abundance is en hanced more at higher 
shock velocity (jSchilke et al.lll997t) . It is believed that SiO 
abundance is enhanced as a consequence of grain sputter- 
ing or grain-grain collisions releasing Si-bearing material 
into the gas phase, which reacts rapi dly with O-bearing 
species (e.g., Q 2 an d OH) to form SiO (jSchilke et al.lll997t 
iCaselli et al]ll99l . The emission is consistent with its 
production in C-type shocks, with a total velocity range 
of ~ 20 km s -1 (see Fig. [T0|) . similar to that found to 
produce the observed SiO c olumn densities in molecular 
outflows ( Sc hilke et al.lll997t ). The C-type shocks have also 
bee n argued to produce the high shock velocity gas seen in 
H 2 (jZinnecker. McCaughrean. fc Ravnerlll998h . The SiO 
emission, with the derived density orders of magnitude 
lower than the critical density, may arise from "thinner" 
regions, as in thin shock fronts. 

SiO may trace different shock conditions from H2. As 
mentioned, although H2 emission is bright around bow 
shocks NK1 and SKI, SiO emission is faint around them, 
and is faintest around the former. In H2, the FWHM ve- 
locity widths toward the bow shocks were found to be 9 
(NK7), 6 (NK4), 12 (NK2), 31 (NK1), 21 (SKI), and 6 



(SK2) km s" 1 (|Davis et alj|200d iTakami et alj|2006h . in- 
dicating that the shock is much stronger for bow shocks 
NK1 and SKI, and is strongest for the former. Thus, 
it is possible that SiO around bow shocks NK1 and SKI 
is mostly destroyed by the shock s through reactions such 
as SiO+OH^Si0 2 +H (see, e.g JSchilke et al.lll997ft . On 
the other hand, SiO and H2 are known to trace different 
physical conditions, with SiO (Tex ~ 100 K, rtcr ~ 10 8 
cm -3 ) tracing much cooler but much denser gas than H2 
(T ex ~ 2000 K, n cr ~ 10 6 cm -3 for H 2 -H 2 collision). 
Thus, it is also possible that although gas-phase SiO is be- 
ing produced, the density there is not high enough to pro- 
duce bright collisionally excited SiO emission. CO emis- 
sion (Tex ~ 50 K, n C r ~ 4x 10 4 cm -3 ), which traces cooler 
and less dense gas than SiO emission, is seen around there, 
also supporting this possibility. This possibility may ap- 
ply as well to weak bow shocks such as NK6, NK7, and 
SK5. For those weak bow shocks, however, it is also pos- 
sible that SiO there is depleted back onto grains in the 
postshock gas. 

The density in the bow shocks is expected to decrease 
with the distance from the source due to the sideways ejec- 
tion of the shocked material and the possible thermal ex- 
pansion of the jet. The fact that the CO emission de- 
creases with the distance from the source (see Fig. ITOl) 
also supports this possibility. However, due to the shock 
enhancement, the SiO emission is brightest around bow 
shocks NK4 and SK2. 

CO emission is also seen tracing the shocks in the jet. 
For those bow shocks, e.g., NK4 and SK4, where CO 
and SiO emission are coincident, the velocity dispersion is 
smaller in CO than SiO, suggesting that CO traces weaker 
shock than SiO. It is likely because CO traces cooler and 
less dense gas than SiO, and thus further away from the 
bow tip where the shock velocity is lower. 

4.1.2. Inclination 

The jet inclination, i, can be estimated from the mean 
velocity of the SiO emission, v m , with i = sin" (v m /vj), 
where Vj is the jet velocity. The northern component and 
southern component of the jet, however, may have differ- 
ent inclinations because of their different mean velocities. 
Assuming the same jet velocity, the northern component 
and southern component of the jet have an inclination of 
- 6.3°(100 km s~ X /vj) and 4.0°(100 km s^/vj), respec- 
tively, in agreement w ith that found from the water masers 
(jClaussen et al.lll998l ). Note that, knots SS and SN may 
have lower inclinations, because of their lower mean veloc- 
ities. 

4.1.3. Mass-loss rate 

Continuous structures between the bow shocks may 
trace the "intrinsic" jet itself, allowing us to estimate the 
mass-loss rate of the jet, which is given by 



Mi 



TTd 2 



-Vjmu 2 



(1) 



where n c and d c are the number density and the diameter 
of the continuous structures. With n c ~ 5.8 x 10 4 cm -3 , 
d c ~ 460 AU, (i.e., 1"), and v 3 = 100 - 200 km s -1 , the 
mass-loss rate is (1 - 2) x 10~ 6 M Q yr -1 , ab out 15%- 30% 
of the infall rate derived from the envelope ([Paper J ). 
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4.1.4. Episodic nature 

In the inner part of the jet, a chain of bow shocks are 
seen with a semiperiodic spacing of ~ 3". Between the bow 
shocks, the continuous structures are seen probably with 
the velocity increasing periodically from one bow shock 
to the next, suggesting that the jet velocity varies periodi- 
cally with time and is highest at the far ends of the contin- 
uous structures. These are consistent with a pulsed (i.e., 
episodic) jet model with a periodical velocity variation, in 
which a chain of bow shocks are formed as the fast-moving 
jet material impacts on t he slow-moving jet material (see, 
e.g.. ISuttner et al.lll997t iLee et alj|200ll ). The jet like SO 
emission around the source, where the velocity probably 
has not steeped into shock, may show intrinsically how the 
velocity varies with time. Observations at higher angular 
resolution are needed to resolve it. 

The period of the velocity variation, which can be es- 
timated by dividing the semiperiodic spacing between 
the bow shocks by the jet velocity, is found to be P ~ 
65(100 km s^ 1 /vj) yr. The jet is believed to be launched 
from an accretion disk around the source. The periodic 
velocity variation may be due to periodic perturbation of 
the accretion disk in an eccentric binary system; in this 
case the period of the velocity variation is equal to the or- 
bital period of the binary and the separation of the binary 
is given by 

a=(^r (2) 

wher e M is th e total mass of the binary. With M ~ 0.15 
Mq ([Paper J ) and P = 65 yr, we have a ~ 8 AU (i.e., 
0^02) , similar to that estimated in HH 34 (jReipurth et alJ 
2002). A binary with this separation, however, would be 
unresolved in our observations. The periodic variation 
may also be due to FU Ori like events as the disk builds up 
mass from the envelope then dumps it onto the protostar 
in a burst. 

4.1.5. Precessing 

The jet itself may be slightly precessing as the jetlike 
SO emission shows a slight S-shaped symmetry about the 
source. The fact that the continuous structures seen in H2 
and CO are sinuous also supports this possibility. The an- 
gle of precession can be assumed to be the angle subtended 
by the continuous structures and is found to be < 2°. Since 
the continuous structures may show a full cycle between 
bow shocks SKI and SK2, and between bow shocks NK1 
and NK2, the period of the precession is probably similar 
to that of the velocity variation, suggesting that the pre- 
cession is also due to a binary system. A slight S-shaped 
symmetry is also seen in other jets, e.g., HH 34, and is 
generally ascribed to the tidal effe cts of the companio n 
star on the direction of the jet axis (jReipurth et a l. 2002). 
However, it is also possible that the slight preces sion is 
due to kink instability in the jet (jTodo et al.l 19931) 



4.1.6. Bending 

The inner part of the jet is bent by ~ 1° to the 
west. Similar jet bend ing is also seen in HH 211 
(jGueth fc Guilloteaul [T999). It has been proposed that a 
jet can be bent (or de flected) by the dynamical pressure of 
the ambient medium (jFendt fc Zinneckerlll998h . However, 



it is not clear how the ambient medium can be communi- 
cated to the highly supersonic jet. 

The bending c ould be due to motion of the jet source in 
a binary system (jFendt fc Zinneckerl 1998). In this case, a 
binary system with a minimum separation of ~ 0"5 (or 230 
AU) , which is the lower limit set by the half displacement 
of bow shocks NB1/2 and SB 1/2 from the original axis, 
is required. However, no binary companion has been de- 
tected yet. It is possible th at the bending is due to L orentz 
force on the magnetic jet (jFendt fc Zinneckenll998l ). It is 
also possible that the jet can appear to be bent just due 
to the motion of the star/disk system, i.e. the jet is orig- 
inating from different positions as a function of time due 
to the star's peculiar velocity. 

4.1.7. Rotation? 

In theoretical jet-launching models, the jet is expected 
to be rotating, carrying away angular momentum from the 
infalling envelope and the accretion disk. The far end of 
the southern component of the SO jetlike structure, which 
is at a beam size (~ 1") away from the source, is likely 
far enough from the possible disk contamination around 
the source and thus can be used to study the jet rota- 
tion. However, no clear velocity gradient is seen across 
it (Fig. [TTa). Knots SS and SN are the shocks closest 
to the source and can also be used to study the jet rota- 
tion. Although no clear velocity gradient is seen across 
knot SS (Fig. [TTb). a velocity gradient is seen across knot 
SN (Fig. [TTb ) with the redshifted emission to the east and 
the blueshifted emission to the west, similar to that seen 
in the envelope. Thus, jet rotation may be seen in our ob- 
servations. Observations at higher angular resolution are 
needed to confirm this. 

4.2. Molecular Outflow Shells 

Outflow shells are seen in CO J = 3 — 2 surrounding 
the jet extending to the north and south from the diffuse 
nebulae that trace the outflow cavity walls (M ZAML02), 
similar to that seen in CO J = 2 - 1 dPaper J ). However, 
the shells are asymmetric in CO J = 3 — 2, with the north- 
ern shell bright in the west and southern shell bright in the 
east. Since CO J = 3 — 2 traces warmer material and thus 
more recent outflow interactions than CO J = 2 — 1, this 
asymmetry indicates that the recent outflow interactions 
are asymmetric. This can be explained if the jet axis has 
rotated in the clockwise direction recently, as suspected 
from the morphology of the H2 jet (see £|3.2p . 

Around the source, internal outflow shells are seen in 
CO and HCO + associated with, e.g., bow shocks SKI 
and NK1. They probably trace the material squirted out 
from the bow shoc ks, as in pulsed jet model (see, e.g., 
ISuttner et aHll997HLee et al.ll200lh . The shells are asym- 
metric, with the northern shell bright in the east and the 
southern shell bright in the west. This is expected if the 
shells are produced by the slight S-shaped jetlike SO struc- 
ture. The bases of the HCO + shells, which are seen with 
a hint of rotation, however, are probably consisted mainly 
of the material extended from the inner envelope and even 
the possible disk. The bases are also seen in H 13 CO + and 
continuum, probably tracing the dense material extended 
from around the same regions. 

These internal outflow bases may trace the material en- 
trained from the inner envelope by the jet. However, in the 
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currently existing jet simulations (see, e.g.. ISuttner et al 
1997t iDownes k Ravi Il999t iLee et alj 120011: iRaga et al 



2004), no internal outflow bases have been seen because 
the jet is protected by a slower sheath shock from con- 
tinued interaction with the ambient material at its bases. 
It is possible that these internal outflow bases trace the 
material directly launched from the inner envelope (and 
even the possible disk) around the jet and then swept up 
by the internal bow shocks, e.g., NK1 and SKI, of the jet. 
Further simulations are needed to study this. 

4.3. Wide-angle wind component? 

The jet itself is clearly highly collimated, yet the diffuse 
nebulae clearly show the classic bipolar parabolic shape 
with a wide opening angle. This may suggest that an un- 
seen wide-angle wind component is there to produce the 
diffuse nebulae, so that the jet is only a dense component 
of a wider wind. However, as seen in Figure [2^,, the jet 
produces a chain of bow shocks with the interacting sur- 
face growing with the distance from the source, and thus 
producing a few prominent bow shocks with much bigger 
interacting surfaces than the jet itself. Through these bow 
shocks, a jet may be able to produce wide-opening outflow 
bases such as the diffuse nebulae. Detailed modelings are 
needed to check this. 

4.4. Molecular Envelope and Disk? 

The HCO + emission around the source may arise from 
the inner part of a dynamically infalling envelope with ro- 
tation, showing a kinematics similar to that seen in the 
C 18 and 13 CO em ission in iPaper J . but closer to the 
source. As argued in iPaperll . a compact rotationally sup- 
ported disk with a radius of ~ 74 AU (or 0"16) is expected 
to be formed within the inner part of the envelope. In or- 
der to study it, the PV diagrams of HCO+, CO, and SO 
emission cut across the source are compared with a Keple - 
rian rotation law (with a source mass of 0.15 M^. lPaper J ) 
and a rotation law with a constant specific angular momen- 
tum (Fig. [T2"|) . In the diagrams, the large-scale emission 
structures are excluded (see in order to study the in- 
ner regions. Note that the comparison could be affected 
by the missing flux around the systemic velocity. Rotation 
is clearly seen within 100 AU from the source in CO and 
HCO + , with the redshifted emission on the east and the 
blucshifted emission on the west. The SO emission may 
have two components, with the low velocity (< 2 km s _1 
from the systemic velocity) from the rotation and the high 



velocity from the jet. The PV structures seem better fit- 
ted by the rotation law with a constant specific angular 
momentum, indicating that the emission within 100 AU 
from the source probably is still from the inner part of a 
dynamically infalling envelope with rotation rather than a 
rotationally supported disk. The rotation may change to 
Keplerian in the innermost region. Observations at higher 
angular resolution and detailed modeling are both needed 
to check this. 

5. CONCLUSIONS 

We have mapped the protostellar jet HH 212 in 850 /im 
continuum, SiO J = 8-7, CO J = 3-2, SO Nj = 8 9 -7 8 , 
HCO+ J = 4-3, and H 13 CO+ J = 4 - 3 emission. Ther- 
mal dust emission is seen in continuum around the source 
IRAS 05413-0104, mainly arising from an inner envelope 
and a possible disk. The inner envelope is also seen with 
rotation in CO, HCO + , and probably SO. The structure 
is unresolved but likely to be flattened perpendicular to 
the jet axis, as seen in the CLEAN component map. Like 
H2 emission, CO and SiO emission are seen along the jet 
axis but extending closer to the source, tracing the bow 
shocks with a broad range of velocities and the continuous 
structures in between. SO emission is seen only around 
the source, forming a jetlike structure extending along the 
jet axis from the source, likely tracing the jet near the 
launching region. The jet is episodic and bending. It may 
also be slightly precessing as the jetlike SO structure shows 
a slight S-shaped symmetry about the source. A hint of 
jet rotation is also seen across the jet axis. Internal out- 
flow shells are seen in CO and HCO + , associated with the 
bow shocks in the inner part of the jet. The bases of the 
HCO + shells are seen with a hint of rotation similar to 
that seen in the inner envelope, probably consisted mainly 
of the material extended from the inner envelope and even 
the possible disk. The bases of the outflow shells arc also 
seen in H 13 CO + and even the continuum, probably trac- 
ing the dense material extended from around the same 
regions. Outflow shells are also seen in CO surrounding 
the jet extending out from the H2 nebulae. They are prob- 
ably internal outflow shells associated with the bow shocks 
further down the jet axis. 
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Table 1 



Line properties 



Line 


Frequency 


Eup 


n cr (50-100 K) 




(GHz) ' 


(K) 


(cm -3 ) 


SiO J = 8 - 7 


347.330631 


75.0 


1.1x10 s 


SO Nj = 8 9 - 7 8 


346.528481 


78.8 


1.4xl0 7 


CO J = 3-2 


345.795991 


33.2 


3.6xl0 4 


HCO+ J = 4-3 


356.734288 


42.8 


9.0xl0 6 


H 13 CO+J = 4-3 


346.998338 


41.6 


8.1xl0 6 



n cr : Critical density of molecular hydrogen. 
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Offset (") v (Hz) 

Fig. 1. — (a) 850 (im continuum contours plotted on top of the H2 image adopted from MZAML02. Note that the image has been 
rotated by 22.5° clockwise. Gray contours are from the restored map with a synthesized beam of l"16x0"84, going from 4 to 40 cr with a 
step of 4 cr, where a = 5 mjy beam -1 . Black contours are from the CLEAN component map. The cross marks the source position, (b) 
A single-temperature fit to the spect ral ener gy distribution of the continuum source (see text for detail). The squares are from our SMA 
observa tions, with the open o ne from lPaper II and the filled one from this paper. The open circles, triangle, and diamond are from the IRAS, 
JCMT l|Zinnecker et al.lfl992T ). and VLA l|Galvan-Madrid et aI.ll2004T ) observations, respectively. 



Fig. 2. — (a) The H2 image adopted from MZAML02. The ellipse outlines the observed region in our observations. The dotted lines 
indicate the jet axes. The cross marks the source position, (b) The H2 image rotated by 22.5° clockwise, (c) SiO (integrated from -21.1 to 
16.0 km s — 1 ) contours plotted on top of the H2 image. The contours go from 3 to 21 cr with a step of 2 cr, where a = 0.58 Jy beam -1 km s -1 . 
(d) High-velocity CO (integrated from -18.4 to -7.1 and from 3.4 to 13.3 km s -1 ) contours plotted on top of the H2 image. The contours go 
from 3 to 21 cr with a step of 3 cr, where cr = 0.38 Jy beam -1 km s . The synthesized beams are 0"96x0"69 and l"16x0"84, respectively, 
for the SiO and CO emission. 
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2 -2 2 -2 2 -2 

Offset (") 

Fig. 3. — (a) SiO, (6) High-velocity CO, and (c) SO contours on top of the H2 image near the source. The SiO and CO contour levels are 
the same as those in Fig. and[2Jl, respectively. The SO (integrated from -11.8 to 11.1 km s~ 1 ) contours go from 2 to 10 a with a step of 
2 cr, where cr = 0.45 Jy beam -1 km s -1 . The synthesized beams are 0!'96x0!'69 for the SiO and SO emission and l"16x0"84 for the CO 
emission. 
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10 -10 10 -10 4 2 -2 -4 

Offset (") 

Fig. 4. — CO contours plotted on top of the H2 image. The synthesized beams all have a size of l"16x0"84. (a) shows the total CO 
emission integrated from -18.4 to 13.3 km s _1 . The contours go from 4 to 36 a with a step of 4 cr, where cr = 0.55 Jy beam -1 km s _1 . (6) 
and (c) show the low-velocity CO emission integrated from —4.0 to 3.1 km s _1 . The contours go from 5 to 45 cr with a step of 5 cr, where 
cr = 0.24 Jy beam -1 km s _1 . 
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Velocity offset from systemic (krn/s) 

Fig. 5. — PV diagrams of the SiO and CO emission cut along the jet axis with a width of 0"9. (a) SiO (gray) and high-velocity CO (green) 
contours on top of the H2 image. (6) PV diagrams of the SiO (black contours with image) and CO (green contours) emission. Here, HV and 
LV denote the high-velocity and low-velocity ranges, respectively, for the CO emission. The red and blue dashed lines mark the near and far 
ends of the continuous structures. 
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Fig. 6. — PV diagrams of the SiO, CO, and SO emission cut along the jet axis with a width of 0"9. (a) SiO (black) and SO (red) contours. 
(b) PV diagrams of the SiO (black contours with image), CO (green contours), and SO (red contours) emission. 
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2 -2 2 -2 2 -2 2 -2 

Offset (") 

Fig. 7. — HCO + , continuum, and H 13 CO + contours plotted on top of the H2 image. The synthesized beams all have a size of l"16x0"84. 
(a) HCO -1- (black, integrated from -2.5 to 2.5 km s — 1 ) and low-velocity CO (green) contours. The HCO -1- contours go from 3 to 15 a with a 
step of 2 cr, where a = 0.32 Jy beam -1 km s — . (b) Blueshifted (integrated from -2.5 to km s —1 ) and redshifted (integrated from to 2.5 
km s — 1 ) HCO + emission, with the contours from 3 to 11 and 3 to 9 cr, respectively, with a step of 2 <r, where a = 0.23 Jy beam -1 km s — . 
(c) HCO + (black) and continuum (orange) contours, (d) HCO + (black) and H 13 CO" 1 " (cyan, integrated from -1.7 to 2.5 km s — 1 ) contours. 
The H 13 CO + contours go from 3 to 6 a with a step of 1 cr, where a = 0.20 Jy beam -1 km s — 1 . The horizontal dashed lines indicate the cuts 
for the PV diagrams in Fig. [8] 
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Fig. 8. — PV diagrams of the HCO+ (black) and H 13 CO + (gray) emission cut across the jet axis, with the cuts shown in Fig[7}i. (a) A cut 
centered at the source. Cuts centered at (f>) 2"5 to the north, (c) 1"25 to the north, and (d) 1"25 to the south, respectively, of the source. 
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Fig. 9. — (a) Line ratio of CO J = 3 — 2/ J = 2 — 1 smoothed to the angular resolution in the CO J = 2 — 1 observations, which is ~ 2"5. 
(b) Excitation temperature derived from the line ratio. The black contours in these two panels show the PV diagram of the CO J = 3 — 2 
emission cut along the jet axis, smoothed to the resolution of the CO J = 2 — 1 observations. 
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Fig. 11. — PV diagrams with cuts across (a) the southern component of the SO jetlike structure, (6) knot SS, and (c) knot SN. The cuts 
have a width of 0"15. The dashed vertical lines indicate the possible jet central velocities for the cuts. In (6) and (c), the black contours are 
from the SiO emission and the gray-scale image with gray contours is from the CO emission. 
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Fig. 12. — PV diagrams of CO (green), HCO + (red), and SO (blue) emission cut across the source. The dotted curves are calculated 
assuming a Keplerian rotation with a source mass of 0.15 Mq. The dashed curves are calculated assuming a constant specific angular 
momentum with v ro ^ = 2^p- k m s -i j with d being the projected distance in arcsec. Note that the fits could be affected by the missing flux 
around the systemic velocity. 
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